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Abstract 

We present results on the stability of quantum systems consisting of a nega- 
tive charge —qi with mass mi and two positive charges q^ and (73, with masses 
m 2 and 777,3, respectively. We show that, for given masses mj, each instability 
domain is convex in the plane of the variables {qi/q2, Qi/Qz)- A new proof 
is given of the instability of muonic ions (a,p,fi~). We then study stability 
in some critical regimes where qs <C qi'- stability is sometimes restricted to 
large values of some mass ratios; the behaviour of the stability frontier is 
established to leading order in q^/q2- Finally we present some conjectures 
about the shape of the stability domain, both for given masses and varying 
charges, and for given charges and varying masses. 

' Dedicated to the memory of Harry Lehmann 

* Work supported in part by the U.S. Department of Energy under Grant 
No. DE-FG02-84-ER40158. 

CERN-TH 98-351 
February 2, 2008. 



1 Introduction 



In two previous papers (T|, 0, hereafter referred to as I and II, respectively, 
we studied the stability of quantum systems consisting of point-like electric 
charges 

Qi = ±[-91,92,93], a > 0, (1) 

and masses m«. The Hamiltonian is 

H = -S- + + -S- - ^ - + M (2) 

2mi 2m 2 2m 3 r 12 r 13 r 23 ' 

where = and = \ fi — fj\. 

In I, we considered the case of unit charges = 1, with application to 
physical systems such as H 2 + (e~pp), H~(pe~e~) and Ps~(e + e~e~), which 
are stable, or (e~pe + ) which is unbound. We pointed out simple properties 
of the stability domain, leading to a unified presentation of results already 
known ||, and to a number of new results. For instance, if m is the largest 
proton mass which gives a stable (e~,e + ,p) system when associated to a 
positron and an electron both of mass m e = 1, it is found in I that m < 4.2, 
a significant improvement over previous bounds ||. This means that global 
considerations on the stability domain can sometimes complement specific 
studies adapted to particular mass configurations. 

In II, we extended the discussion by letting the charges g, themselves 
vary. The number of parameters is increased from two to four, and one can 
choose two mass ratios and two charge ratios. The general properties of the 
stability domain established in II will be briefly reviewed, and supplemented 
by new results. This will be the subject of Sec. ^|. 

As an example of application of general considerations on the stability 
domain, we shall present in Sec. ^ a new proof that muonic ions involving 
a helium nucleus a, such as (p,a,fj,~) or (d, a, //~), are not stable. This 
confirms results Q obtained previously using the Born-Oppenheimer frame- 
work. 

In II, we also considered the limiting case where g 3 — * 0, but only in the 
Born-Oppenheimer case where m 2 = m 3 = 00. We shall resume in Sec. f| our 
investigations and study the behaviour of the stability frontier in the case 
where q 2 > 1 and g 3 <ti g 2 - 

In Sec. ||, we shall present some speculations about the plausible shape 
of the stability domain in both representations: varying charge-ratios for 
given masses, or varying masses for given charges. A number of interesting 
questions remain open. 

Our rigorous results are supplemented by numerical investigations based 
on a variational approximation to the solution of the 3-body Schrodinger 
equation. In particular, we display an estimate of the domain of stability in 
the (g 2 , g 3 ) plane for some given sets of constituent masses. 
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2 General properties of the stability domain 



2.1 Inverse-mass plane for unit charges 

Consider first the case where q\ = qi = 53, which can be chosen as g» = 1. 
Stability is defined as the existence of a normalised 3-body bound state with 
an energy below that of the lowest (1,2) or (1,3) atom, i.e., 

< min (E?l E®) , E$ = -(a 1 + a,)" 1 ^, (3) 

where = 1/mj is the inverse mass of particle i. Thanks to scaling, there 
are only two independent mass ratios in this problem. In I, we found it 
convenient to represent any possible system as a point inside the triangle of 
inverse masses normalised by Ylt=i a i — ^- This triangular plot is shown in 
Fig. § 

In this representation, the stability domain appears as a band around the 
symmetry axis where «2 = 03. It is schematically pictured in Fig. ffl. 




Figure 1: Domain of possible inverse masses aj, normalised by J^ctj = 1, 
and shape of the stability domain for three unit charges. 

The following rigorous properties are known, or shown in I. 

a) All points of the symmetry axis («2 = 03, < a% < 1) belong to the 
stability domain ||. 

b) Each instability domain is star shaped with respect to the vertex it con- 
tains. For instance, at the right-hand side of Fig. [I], each straight line issued 
from A 2 crosses at most once the stability frontier between A 2 and the sym- 
metry axis. 

c) Each instability domain is convex. 
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2.2 Inverse-mass plane for unequal charges 



For arbitrary charges % the threshold energy (|3]) is modified 

" 2{a 1 + a i y W 

The separation between the two thresholds, (T), which plays a crucial role 
in the discussion, is given by 

Oi + a 3 )q% = Oi + a 2 )ql- (5) 

It is a straight line in both pictures, i.e., for fixed charges qi in the plane of 
inverse masses, and in the charge plane for fixed masses. 

One expects an increase of stability near (T), where both thresholds be- 
come equal. This is what happens in the unit-charge case where, according 
to Hill's theorem ||, we have stability for a 2 = a%. Another example is the 
Born-Oppenheimer limit with m<i = 777,3 — 00, and say ni\ = q\ = 1 to fix 
the scales: in the (q 2 , 93) plane, it is observed that the stability domain does 
not extend much beyond the unit square (q 2 < 1, #3 < 1), except for a spike 
around the q 2 = q% axis, which reaches q 2 — (fe — 1.24 Q. This is shown in 
Fig. 1 




Figure 2: Schematic shape of the stability domain in the Born-Oppenheimer 
limit. The heavy particles have charges q 2 and ^3. The charge of the light 
particle is set to —q\ = — 1. 

For fixed charged qt, the triangular plot of Fig. p] can be used again. The 
threshold separation (T) is a straight line passing through the (unphysical) 
point oti = — 1, a 2 — CU3 = 1, which is the mirror image of Ai with respect to 
A2A3. As seen in II, each instability domain remains star shaped and convex, 
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as in the case of unit charges. In particular, if g 3 < q 2 and q 3 < 1, the entire 
sub-triangle limited by (T) and A 2 corresponds to stability. If furthermore 
q 2 < 1, then we have stability everywhere. 

Following a suggestion by Gribov |7|], one can also consider level lines of 
constant relative binding, i.e., such that 

E[f<E$ and E®/E$=\>1. (6) 

(Remember that both and E^ 2 are negative.) What happens is that 
the domain E^/E^ < A (including the case where particle 3 is unbound, 
where we set 

£0) = e$) is also convex and star shaped. 
The proof is essentially the same as for the domain of instability. One 
first rescales the ctj from a± + a 2 + a 3 = 1 to a.\ + a 2 = 1, so that the 
threshold energy E^ 2 becomes constant. If two points a = (ai,a 2 ,a 3 ) and 
a 1 = (a[, a 2 , a' 3 ) belong to the frontier of interest, that is to say 

E^(a) = E^(a') = XEil (7) 

then as the a, enter the Hamiltonian linearly, for any intermediate point 
xa + (1 — x)a' with < x < 1, one has 

E^\xa + (1 - x)a') < xE (3 \a) + (1 - x)E^\a') = \E$ . (8) 

Similarly, a decrease of a 3 with a,\ and a 2 kept constant cannot do anything 
but decrease E^ with E^ unchanged: this proves the star-shape behaviour. 

2.3 Convexity in the , q% variables 

We return to the domain of strict instability, but now for fixed masses 
and variables charges. We fix qi — 1 and consider the frontier of stability in 
the (l/<?2, 1/93) plane. 

First, we notice that the domain of stability is star-shaped with respect 
to the origin. Indeed, when a system of charges (— l,g 2 ,g 3 ) is transformed 
into (— 1, g 2 /A 2 , g 3 /A 2 ), with A > 1, the new system can be rescaled into 
(— A, g 2 /A, g 3 /A) which experiences the same attraction but less repulsion 
than the original system. 

If m 2 < m 3 , the threshold separation (T) has a slope larger than unity 
in the (I/92, 1/93) plane. If m 2 > m 3 , this is the reverse. Consider now for 
definiteness the domain where 

q\ («! + a 3 ) > ql {oi X + a 2 ) . (9) 

In this domain, the (1,2) atom is more bound than (1,3) and it is the energy 
of (1,2) to which E^ should be compared. The ground-state energy E^ of 
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the Hamiltonian (^) is separately and globally concave in q 12 , q 23 and g 2 3- 
With our choice of the charges we have 



912 = 92, 



913 = 93, 



923 = 9293, 



(10) 



but we can make a rescaling taking 



912 = 1, 




923 = 93- 



(11) 



In this way, the binding energy of (1,2) is fixed. As is a concave function 
of gi 3 and g 2 3, the instability domain is convex in the (913,923) plane. This 
means a segment of straight line 



joining two points of the domain also belongs to the domain. But this equa- 
tion translates into 



and thus also represents a straight line in the (I/92, 1/93) plane. Thus each 
instability domain is convex in this plane. This is schematically pictured in 



3 Instability of \i ) systems 

The problem of the instability of ions such as (a,p,fi~), (a,d,e~), etc., has 
been considered by several authors. What has been shown is essentially that, 
within the Born-Oppenheimer approximation, the effective (a,p) or (a,d) 
potential is unable to support a bound state ||. The proof below is more 
general, for none of the masses is assumed to be very large. 

First we notice that for equal charges q 2 = 53 > 1.24 || the 3-body system 
is unstable for equal masses 0:2 = 0:3, because it is unstable for a 2 = a 3 = 
and a 2 = a 3 = 1/2, and the domain of instability with respect to a given 
threshold is convex in the triangle of inverse masses. 

Furthermore, using for instance the (52, 93) plane and keeping the masses 
equal and constant, we know that the system is unstable for a fixed q 2 > 1.24 
and q 3 very small, because particle 3 is almost free and because the system 
(1,2) is repulsive. So using convexity in q 3 for fixed q 2 , we prove that the 
system is unstable for < q 3 < q 2 , if 92 > 1.24, for a 2 = a 3 . 

Therefore, if a 2 = a 3 , a system with q 2 = 2 and ^3 = 1 is unstable. But 
the system (1,2) is repulsive and for a.\ = a 2 = 0, i.e., mi = m 2 = 00, we 
have instability. Now, in the left half-triangle a 2 < a 3 , where 



9i3 = ^923 + b 



(12) 




(13) 



Fig. B 



(a± + 0:3)^2 = 4(«i + a 3 ) > («i + ^2)93 



(14) 
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Figure 3: Shape of the stability domain in the plane of inverse charges 
(g^ -1 ,^ 1 ), with normalisation q\ = 1. The lines q 2 x = 1 and q 2 x = 1 are 
either asymptotes, or part of the border, starting from a value of q 2 x or q^ 1 
which might be less than 1, unlike the case shown in this figure. 

so that it is the (1, 2) system which is more negatively bound. Therefore, 
we can use the star-shape instability in the whole left half-triangle, which 
includes not only apfi~ or ape~ , but also adfi~ or odjjT . Notice that the 
proof does not work for q 2 > q%, m 3 > m 2 . 

4 Some limiting configurations 

4.1 The very asymmetric Born— Oppenheimer case 

In II the case was considered where the masses m 2 and m 3 are both infinite, 
but the corresponding charges q 2 and q$ are freely varying. The stability 
domain is shown in Fig. ^, with a normalisation qi = 1. The domain is of 
course symmetric under (g 2 <-> q%) exchange, and includes the (q 2 < 1, (ft < 1) 
unit square. We already mentioned the peak at q 2 = q% — 1.24. The frontier 
starts from q 2 = l,q 3 = 0, where a behaviour 



- 1 ~ 



18 



(-ln<? 3 ) 3 ' 



(15) 
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is proved in II. This leading order ( |T5D is however very crude, as the first 
corrections differ only by terms with higher power of (— lng 3 ). 



4.2 The Born— Oppenheimer approximation for mi and 
m 2 very large, q t = q 2 = 1 

We now consider systems analogous to (p,p,e~). In the limit of strictly 
infinite masses, we have a point source of charge q 2 — q\ acting on the charge 
q% > 0. Thus there is no binding if q 2 > q\ and binding for q 2 < q%. So, a 
non-trivial case consists of mi and m 2 very large but finite, and q\ = q 2 = 1. 
We argue below that binding is unlikely if m 3 g 3 is very small. 

With obvious notations, the adiabatic approximation relies on the decom- 
position 

H = --— + h(r) - 1/r 

2/^12 

L Aa (16) 



2m 3 r 13 r 23 

with = fn\ X + fn 2 ■ Then H > H, where H is deduced from H by 
replacing h by its ground-state energy or its infimum, say inf(h), which is a 
function of r = r 12 . 

The very crude inequality 



r 23 < r + n 3 (17) 



leads to 

As it is known that 

we are sure that 
On the other hand, 



h -~^-^- (18) 



A3 - -X- > 0, (19) 

^'13 



h > if 8m 3 g 3 r < 1. (20) 



> -*L - «L > (21) 
- 2m 3 r 13 " 2 1 ; 



which is obtained for r — > oo. So 

f = for r < (Smsgs)" 1 = R 

ml[h) | > _ m3 g2/ 2 for r > i?. 

(Of course, inf h(r) must be a continuous function leaving for some r > R, 
and reaching — m 3 gf/2 at large r.) 
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Consider now the relative motion as described by 

H 



~ + in£(h)--. 
2fx 12 r 



(23) 



If we neglect mi(h), H reduces to the Schrodinger equation for a two-body 
atom, with ground-state energy — /i 12 /2 and reduced radial wave function 
u(r) = 2j2% 2 r exp(—jj,i 2 r) . For /i 12 ^> 1, this wave function is concentrated 
near r = and inf(/i) can be considered as a perturbation. The first order 
correction 5E is negative and is such that 



\SE\ < 



m z ql 



u(r) 2 dr 



A.2 + ^12^3 ra 3 ql 
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H 12 


exp 


4m 3 q 3 _ 



(24) 



R 



In other words, SE vanishes exponentially when mi and m 2 — > oo for fixed 
(J3 and 7723. This is beyond the accuracy of the Born-Oppenheimer approxi- 
mation, and strongly suggests that there is no binding if 7773(73 <C // 12 , which 
reduces, for mi = 00, to m 3 q 3 <C m 2 . 



4.3 The small m 3 q s limit for mi = 00 

We now extend our study of the configurations with charges g 2 > 1 and 
<?3 ^ 92, with normalisation g x = 1. Instead of the Born-Oppenheimer limit, 
we consider the somewhat opposite case where m x = 00, i.e., the lower side 
A 2 A 3 of the triangular plot. The threshold separation (T) happens for 

«3 ( <?3 X 1 



. . (25) 

«2 \<?2/ 

very close to A 2 . Some crude variational calculations has convinced us that 
the frontier occurs for a 3 /a 2 = 0(q 3 ), not 0(q 3 ). This suggests a first order 
calculation in q 3 . We temporarily fix the scale at a 2 = 1 and split the 
Hamiltonian into 

H = H„ + q 3 H u 

H„ = -\^~% (26) 



2 V^3/ ^3 ^23 

where ?i = 0, and r 23 = |r 3 — r 2 |. We are faced with a standard problem of 
degenerate perturbation theory. At zeroth order, we get the energy E and 
eigenfunction ty 

£o = -f, vD = V,(r 2 Mr 3 ), (27) 
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where ^(r 2 ) = 7r _1 / 2 (g 2 ) _3 ^ 2 exp — (^2^2) and <^(r 3 ), yet unspecified, is deter- 
mined by diagonalising the restriction Hi of Hi to the ground-state eigen- 
space of Ho. This reads 

Hiip(r 3 ) = Ehp(t 3 ) 

/(x) = 1 — (1 + x) exp(— 2x). 

For g2 < 1) the potential in (p8|) exhibits an asymptotic Coulomb be- 
haviour which is attractive. Thus 1 supports bound states whatever inverse 
mass 0:3/(73 is involved. We recover the property seen in (II) that for q% < 1 
and q2 < 1, the 3-body system is stable for any choice of the constituent 
masses. 

For q 2 > 1, the potential in ( fZ8D has a repulsive Coulomb tail or decreases 
exponentially. At best, it offers a short-range pocket of attraction to trap 
the charge ^3. The short-range character is governed by the exponential in 
the form factor /, as per Eq. Such a potential supports a bound state 

provided the mass 53/03 is large enough, say ^3/03 > /x c . This is why at the 
frontier a 3 = 0(qs). 

Calculating the critical mass \i c accurately as a function of q<i is a routine 
numerical work. One can for instance integrate the radial equation at zero 
energy and look at whether or not a node occurs in the radial wave-function 
at finite distance. The result is shown in Fig. [|. 




Figure 4: Minimal reduced mass fi c to achieve binding in the potential 
V — — 1/r + q/r(l — (1 + qr) exp(— 2qr)). 

The behaviour observed in Fig. |4] is not surprising. After rescaling, the 
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Hamiltonian of Eq. (|2"8|) can be rewritten as 
H x = 2q 2 

where q = q 2 and p = q^/a^. The critical mass for achieving binding in 
a Yukawa potential V\ = — exp(— r)/r is well known || and well studied 
||. It is p\ ~ 1.679. For an exponential, it is about 1.446 ||, and thus 
p 2 — 2.892 for V 2 = — exp(— r)/2. It is easily seen that the critical coupling 
p c for binding in V\ + V 2 is such p~ l < p^ 1 + p^ 1 ■ This means p c > 1.06 
for the the attractive part in (P9"D, a bound not very far from the computed 
value /i c ~ 1.10. This corresponds to the case q = 1 in Fig. ^. For g > 1, the 
repulsive Coulomb tail makes it necessary to use a larger value of p c , this 
explaining the rise observed in Fig. |4j. 



1 d 2 
p dr 2 



exp( 



-r) exp( 



1/q 



(29) 



4.4 Stability frontier for small qz/q2, and q 2 > 1 

We just established that for mi = oo, small 93, and 52 > 1, the stability 
frontier lies at some 0:3 ~ qz/p C) where Pcfe) is computable from a simple 
radial equation. We now study how the frontier behaves as mi becomes 
finite. We are near A 2 in the triangle, where a 2 — 1, and ol\ and 0:3 are 
small. 

We introduce the Jacobi variables 

-. r a 2?i + «i?2 

p = r 2 -ri, A = r 3 ■ , (30) 

OL\ + a 2 

in terms of which the relative distances are 

ri2 = P, r 23 = A ■ p, r 3i = -A ■ p, (31) 

«1 + «2 «1 + «2 

and the Hamiltonian reads 



- «i + a 2 A P o "3 + — - — 

2 p 2 \ cti + a 2 



^ + ^ , (32) 

I A + axp/{ax + a 2 )\ |A - a 2 p/(ax + a 2 )| 

besides the centre-of-mass motion, which will be now omitted. 
A first rescaling p — > (ax + ot 2 )p/a 2 results into 

2 \ai + a 2 J ax + a 2 p 

-- (a 3 + A\ - + 4^, (33) 

2\ ai + a 2 / \\ + axp/a 2 \ |A - p\ 
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which is the scale transformed of 



H = --a 2 A 2 - -a 3 A 3 1 — — — — r , (34) 

2 2 r 2 r 23 |r 3 + air 2 /a 2 | 

provided the inverse masses in H are proportional to these in H, and the 
strengths in H to these of H. A convenient rule of transformation of masses 
is 



1 / a 2 



a>2 



a 2 + a 3 \ai + a 2 



2 



a 3 = - ( a 3 + aia2 ) , (35) 

"2 + «3 V ai + a 2 J 

since it changes our triangular normalisation = 1 into a 2 + a 3 = 1. 
For the charges, the simultaneous identification 

n «2 

92 = Cq 2 ■ , 

OC\ + OL 2 

93 = ^93, (36) 

9293 = ^9392, 



results into 



92 = 92, 93 = 93 • (37) 

a 2 



The rescaled Hamiltonian (SI]) slightly differs from the Hamiltonian (|26|) 
corresponding to mi = 00. However, the difference between |r 3 + «ir 2 /a 2 | 
and r 3 is of first order in g 3 , and thus enters at order q 2 in H. We are then 
allowed to write the frontier condition as in the previous section, namely 

« 3 ^ -, (38) 

which, when translated into the original variables, reads, at first order 

cm + "3 — — , (39) 

to be compared with the threshold separation oli + a 3 = g|. This means the 
frontier is at first approximation a straight line, parallel to the side A 3 Ax 
of the triangle of inverse masses, as schematically pictured in Fig. |5[ Note 
that the actual frontier is certainly curved, since the instability domains are 
convex, as reminded in Sec. 0. 
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Figure 5: Expected behaviour of the frontier (F) in a situation where q 2 > 1 
and the ratio q 3 /q 2 is very small, so that the threshold separation (T) is very 
close to A 2 . 

4.5 The small qr 3 limit in the (q 2 , ^3) plane 

Let us consider now the (52, Q3) plane (with q\ = 1) for fixed masses. The 
shape of the stability domain is shown in Fig. || The frontier of stability 
leaves the unit square at some finite value of q 3 . Consider, indeed, q 2 = 1, 
with «i = for simplicity, and a mass scale fixed at a 2 = 1. In a (variational) 
approximation of a (1,2) atom times a function describing the motion of the 

as 



third particle, we can read the calculation of Sec. fO 

1 > q 3 > /x c (l)a 3 ~ 1.10a 3 



(40) 



being a sufficient condition for stability. 

A necessary condition of the same type, i.e., q 3 > fl c (^3 can be obtained 



using the method of Glaser et al. ||10|| . The decomposition 



yields the operator inequality [ |TT|j 
H > H' — [ H 



T23 



T23 



<?2<?3 



r 2 3 



+ q 2 q 3 P(Pr 23 P)- L P, 



(41) 



(42) 



where P is the projector over the ground-state of H = —A 2 — ^2/^2 (times 
the identity in the variable r 3 ). Now H' is the sum of Hq in the variable F2, 
and 



H' 



-a 3 A 3 - q 3 /r 3 + 



(M3 

r 3 



1 - 



x 



-\ -1 



1 + a;/2)exp(-2a;) 



(43) 
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Figure 6: Shape of the stability domain in the (q 2 , q 3 ) plane. The frontier 
of stability leaves the vertical side q2 = 1 of the unit square at some finite 
value of q%. It also leaves the horizontal line q% = 1 at some finite value of 
q 2 , which can either smaller (as in this figure) or larger than 1. 



where x = q 2 r 3l m the variable r 3 . For q 2 = 1, this potential supports 
a bound-state provided g 3 > jl c a 3 , with /2 C > 0.34 from the Jost-Pais- 
Bargmann rule |j~2| , and ji c ~ 0.64 from a numerical calculation (looking 



for nodes in the radial wave function at zero energy). 

If q;i > 0, a reasoning similar to that of Subsec. [O] shows that the 
sufficient condition (HOT) is replaced by 



q 3 > /i c (l)(«i + a 3 ) ~ 1.10(«i + a 3 ), (44) 

where the normalisation is «i + a 2 + 0:3 = 1. 

The result ( PD is of course expected to be better if the computed g 3 is 
small, i.e., if a,\ + a 3 1. 



4.6 Frontier in the (q2,Q3) plane at small 0:3/0:2 

We remain in the (q 2 , q 3 ) plane for fixed masses. We assume a± — for 
simplicity, but some results do not depend on this assumption. We can 
normalise to ol 2 — 1. In the limit where a 3 is small, the threshold separation 
(T), as given by Eq. (p5|), has a very small slope with respect to the q 3 axis. 

The frontier exits out of the unit square at q 2 = 1 and a finite value of 
g 3 which is close to a 3 /i c (l), where /i c (l) — 1.10, according to our previous 
computation. If we look at the frontier outside the unit square, we have two 
questions: 

i) in the lower part of the plot, is the frontier strictly below (T) ? 
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ii) in the upper part, does the frontier overcome the line q 3 = 1 ? 



4.6.1 Lower part of the frontier 

To answer the first question, let us consider a situation where g 3 is close to 

\/2 

but smaller than 1, and thus g 2 ~ « 3 is large. If particles 2 and 3 would 
ignore each other, they would bind around particle 1 with approximately 
the same energy, since we are close to (T), but with different Bohr radii 
namely R3/R2 ~ ct^ 2 <C 1. This suggests the approximation of a localised 
(1,3) source attracting the charge 52, corresponding to a 3-body energy 

* = ql{l - q.)\ (45) 

whose equality with the threshold E% = —q\ gives the approximate frontier 

ll = (46) 
«3 2 - q 3 

which just touches (T) at q 3 = 1, as seen in Fig. [7]. Now this approximation 



93 




Figure 7: Upper bound ( |Hf ) on the lower part of the stability frontier, touch- 
ing the threshold separation (T) for q 3 — 1. A value a 3 /o; 2 = 1/25 is assumed 
here. Note that this bound is not expected to be a good approximation for 
small q2, as it does not delimit a convex domain of instability. 

corresponds to write a decomposition 



H 



H23 + V23 

«3^ 2 93 . 1^2 

T P3 - ^ + 2 P2 



g2(l ~ ^ 
r 2 



+ 92?3 



^23 



1 

r 2 



(47) 
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Figure 8: Schematic shape of the stability domain in the ((72,(73) plane, for 
the limiting case where 03/0:2 is very small. Our constraints cannot exclude 
a very tiny peak along (T) above q 3 = 1. 



and neglect the second term, V23. The spherically-symmetric ground-state 
\l/o of H23 can be chosen as a trial variational wave-function for H. The 
Gauss theorem implies that (^ol^l^o) < 0. Hence the ground state of H 
lies below that of H23, and the actual frontier is below the approximation 
flUI), therefore below (T) as long as (73 < 1. 



4.6.2 Upper part 

We now turn to the question of possible binding above the line q 3 = 1. We 
restrict ourselves to mi = 00, although we suspect that our results are more 
general. Numerical investigations using the method described in Appendix B 
suggest the following pattern. For 777-2 = 1713, a spike is observed on the 
diagonal. It reaches about (72 = Q3 = 1-24 in Fig. |2], corresponding to m-i = 
m,3 = 00, and about q<i = (73 = 1.098 [[L3|] for m 2 — rn 3 C mj. The spike 
remains for moderate values of the mass ratio 7723/777-2, as schematically shown 
in Fig. |^. When, however, 777,3/777,2 exceeds a value which is about 1.8, no spike 
is seen within the accuracy of our calculations, i.e., the frontier seemingly 
coincides with the line (73 = 1, until it reaches (T), as pictured in Fig. [S[ 

We are able to show rigorously below that, for large values of 7773/7772 
no binding occurs above (73 = 1 for (72 < (3/4) (7773/7772 ) 1/ ' 2 . Nothing can 
be said however from this latter value to q<i = (7773/7772) 1//2 on the threshold 
separation (T). In other words, a very tiny peak along (T) overcoming (73 = 1 
cannot be excluded. 

For 777! = 00 and q 3 = 1, the Hamiltonian reduces to 



H 



2777a 



1 

r 3 



+ 



27772 



T2 T 2 3 
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2m 2 r 2 r 23 

Let P be the projector on $(r 3 ), the ground state of H = p|/(2m 3 ) — l/r 3 , 
and \l/ the ground state of H. We have the inequality 

2 



df 3 * — M> 

7" 23 



dr 3 $r 23 $ > 



dr^$ 



(P¥(? a ))' 



(49) 



To estimate ( 1 J r |iif| 1 I r ), we first need (^iPol^)? where each \1/ can be read as 
+ (1 - P)#. We have 

(*P\H \(1-P)V) = 0, 
(^P\H \P^) = (VP\py) (-m 3 /2), (50) 
(tf (1 - P)\H \(1 - P)V) > (tf (1 - P)|(l - P)*) (-m 3 /8), 
and similarly, for ho = p 2 /(2m 2 ) — q 2 /r 2 , 

(*P|/i |(l-P)*) = 0, (51) 
- P)|/i |(l - i 3 )*) > - P)|(l - P)*> (-? 2 2 ^ 2 /2). 

Thus 



where 



(¥|#|tt)> ||P*| 2 (-^) 
+ ||(1-P)^|| 2 

= P 2 /(2m 2 ) - g 2 /r 2 + 



(*P|/i |P*) 



m 3 g 2 m 2 



92 



/ <3> 2 r 23 dr 2 

In a situation where /i does not support any bound state, 

m 3 q\m 2 



[V\H\9) > \\py\\ 2 (-^r) + 1(1 - P)*ll 



and hence 



(V\H\V) > inf 



-m 3 /2 

-m 3 /8 - m 2 q 2 /2. 



(52) 
(53) 

(54) 

(55) 
(56) 



Now the hamiltonian /iq has been studied in Ref . [|10j , and shown not to bind 

*™ < 1.2706. (57) 



Therefore if 



m 3 



m 3 m 3 m 2 g2 
IT < 8 2~ 



and 



2m 2 q 2 
m 3 



< 1.2706, 



(58) 



the system is unstable (the first condition implies — m 3 /2 < — m 2 g|/2, i.e., 
(1,3) is the lowest threshold). For large m 3 /m 2 , the first condition is more 
constraining, so we have no stability above g 3 = 1 from q 2 = to q 2 = 
(3m 3 /(4m 2 )) 1 /2. 
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4.7 Numerical results 

We now display an estimate of the domain of stability in the (q2, qz) plane, 
with normalisation q\ — 1. The method, described in Appendix B, is vari- 
ational. Therefore, the approximate domain drawn here is included in the 
true domain. 

Our investigations correspond to mi = 00, and the mass ratio m 3 /m2 
having the values 1, 1.1, 1.5 and 2. In each case, we show the whole domain, 
and an enlargement of its most interesting part, the spike above q2 = 1 and 
<?3 = 1- 




Figure 9: Variational estimate of the domain of stability for mi = 00 and 
7712/7773 = 1, full view (left) and enlargement of the spike (right). The dotted 
line is the threshold separation (T). 

For rri2 = 777.3, in Fig. || we have a symmetric spike. The location of 
the peak at, q<i = q% = 1.098 reproduces fairly well the values given in the 
literature |13[ . 



For 7712/7713 = 1.1, in Fig. [TU], the spike leaves the horizontal line q% = 1 be- 
fore q2 = 1, as a generalized H~ ion with unit charges and masses (00, 1.1, 1) 



is bound. This is no longer the case for 777,2/7713 = 1.5, as seen in Fig. ITT 



For 7772/7773 = 2, shown in Fig. [12], the domain found in our variational 



calculation is flat. In the lower part, it extends appreciably further than 
indicated by the crude approximation of Eq. ((461). 
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Figure 10: Same as Fig. || for m 2 /m 3 = 1.1 



5 Outlook 



Many questions remain open concerning the stability of 3-charge systems. 
Along the paper, we pointed out that in some limiting cases, more accurate 
results would be desirable. For instance, a question is whether very large 
values of the mass ratio m 3 /m 2 exclude the possibility of binding with q 3 > 1. 

There are also more general questions, concerning domains of some pa- 
rameters for which stability will never be reached, whatever value is given to 
the other parameters. 

For given masses mj, the answer is immediate: there is always a set of 
charges, for instance g 3 < g 2 < q\ = 1, that makes the system stable. 

For given charges g 2 and g 3 , and q\ = 1 to fix the scale, the situation is 
different: one has clearly three possibilities. Region {1} is the unit square 
{q 2 < 1, (fe < 1}, where any mass configuration corresponds to a stable 
ion. Region {2} includes for instance the point g 2 = 2 and g 3 = 0.8: there 
is sometimes stability, mi = oo, m 2 <C m 3 is an example, and sometimes 
breaking into an atom and a charge, as for m\ = oo, m 2 ^> m 3 . Region 
{3} includes points like g 2 = g 3 = 2 for which stability will never occur. 
Determining the properties of the boundary between regions {2} and {3} 
would be very interesting. 

A possible starting point is the result by Lieb [14| , that for a fixed nucleus 
aii = 0, qi — 1, a bound state will never occur if 



1 1 
- + - < 1. 

<72 <?3 



(59) 
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Figure 11: Same as Fig. || for m 2 /m 3 = 1.5 



A simple proof is given in Appendix A. 

This upper bound for possible stability at ra\ = oo (i.e., stability occur- 
ring for at least some value of m 2 /m3) is not too far from the lower bound of 



Fig. [13], obtained from our variational method. More extensive computations 
would be necessary to sketch the shape of the domain of absolute instability, 
in particular by relaxing the condition mi = oo. Note that along the symme- 
try axis, the limit is q 2 = q 3 — 1.098 for mi = oo and m 2 = m 3 finite, while 
it reaches q 2 = q 3 ~ 1.24 for mi finite and m 2 = m 3 = oo. Thus, along the 
symmetry axis, the frontier between regions {2} and {3} is saturated in the 
Born-Oppenheimer limit. On the other hand, for q 2 ^> q 3 or q 2 <C q 3 , the 
question is whether this frontier has q 2 = 1 and q 3 = 1 as actual asymptotes, 



as tentatively pictured in Fig.|TJ] or reached these lines above some values of 
<?2 or q 3 . 
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Appendix A: 

Proof of instability for q^ 1 + q^ 1 < 1 

We give here a proof of the result on instability for all values of m2 and m 3 
if q^ 1 + q^ 1 < 1, provided oi\ = 0, with normalization q\ = 1. 

First, it is shown that rp 2 is a positive operator, in the sense that any 
diagonal matrix element is positive. Indeed, separating the radial and angular 
part of p 2 , 



(^|rp 2 |\I>) 



r*A*d (3) r 
L 2 



^d (3) r- / dQ 



d( 3 )f 



dQ / rdr 



dr 2 

d(rV) 
dr 



rdr 



Now consider the Hamiltonian 
H 



P2 , _Ps_ _ 02 _ 03 Q2qs 



2m 2 2m 3 r 2 r 3 r 23 



(60) 



(61) 
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Figure 13: Variational estimate of the domain of possible stability beyond 
q 2 = 1 and ^3 = 1, for rrii = 00. Inside the domain, stability occurs at least 
for some value of the mass ratio 7712/7713. 



whose thresholds (1,«) with i = 2,3 are governed by the Hamiltonian 

?2 



Pi 



2m,- 



ft 



(62) 



These /ij and the 3-body Hamiltonian if fulfill the identity 

r 2 + r 3 1 



r 3 (H-h 2 )+r 2 (H-h 3 ) = r 2 



_pL 

2m 9 



p1 

-7" 3 - + 92?3 

27773 



1 



^23 



92 ?3 



- " - • (63) 



In the r.h.s., the two first terms are always positive, and so is the third one 
if q 2 x + q^ 1 < 1, due to the triangular inequality. Looking now at the l.h.s., 
its expectation value is always positive, which means that 



[V\r 3 (H - h 2 )\V) > or ($\r 2 (H - h 3 )\^) > 0. 



(64) 



In the first case take \1/ as the ground state of H, which satisfies Hty = E^^. 
This translates into 

£ (3) (^|r 3 |^) > (*|r 3 (■&- - I*), (65) 



2m? 



r-2 



or 



1 / p1 


_ ^ 


V 2777-2 


^2/ 



(66) 
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Figure 14: Guess at the shape of the border of the domain of absolute 
instability. In region {1}, binding occurs for any set of constituent masses. 
In region {2}, binding is achieved under some conditions for the masses. In 
region {3}, stability is never obtained. 



from the variational principle. In the second case 

So, if mi = oo and g 2 <?3 > Q2 + Q3, either E^ > e[^ or E^ > E$ . 



(67) 



Appendix B: Variational method 

We briefly describe the variational method used for the numerical results 
displayed in this paper. More details can be found in . The ground state 
of the Hamiltonian (0) has been searched using trial wave functions of the 



type 16 



^2 C ^ i Pi = y^^Ci [exp( -(^23 - bir 31 - dr 12 ) + 



(68) 



from which all matrix elements can be calculated in close form. The dots 
are meant for similar terms obtained by permutation, in the case of identical 
particles. For given range parameters, the weights C{ are listed in a vector 
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C, which is found, together with the variational energy e, from the matrix 
equation 

f T + V) C = eiVC, (69) 



involving the restrictions of the kinetic and potential energy to the space 
spanned by the tfi, whose scalar products are stored in the positive-definite 
matrix N. 

As the number of terms increases, it quickly becomes impossible to deter- 
mine the best range parameters, even with powerful minimization programs, 



as too many neighboring sets give comparable energies. One way out fl? 
consists of imposing all at, hi and q to be taken in a geometric series. Then 
only the smallest and the largest have to be determined numerically. For 
instance, this method allows one to reproduce the binding energy —0.262005 
of the Ps~ ion, in agreement with the best results in the literature. 

The question now is to find the frontier. Let us consider, for instance, the 
problem of Sec. |4.6.2| . Here mi = oo, q\ = m 2 = 1, and qijqz = m 3 when 



one searches the limit of stability among the threshold separation (T). 
One can estimate the ground-state energy of 

->2 ->2 1/2 1/2 9 

P2 + _Ps_ _ m 3 Qs _ qs + m 3 ' qj 
2 2m 3 r 2 r 3 r 23 

starting from some low value of g 3 , and examine, by suitable interpolation, 
for which g 3 it matches the threshold E t ^ c = —q 3 m^/2. 

A more direct strategy consists to set E = E t ^ e in the Schrodinger equa- 
tion, apply a rescaling and solve 



->2 -2 1/2 i 1/2 

pi pi mJ 1 m 3 mJ 

i^ + Tp + ^r = -<?3 — , 71 

2 2m 3 r 2 r 3 2 r 23 



using the same trial function (|68D , resulting in a matrix equation very similar 
to ( |69|) , where the positive-definite matrix iV now represents the restriction 

1/2 

of m 3 /r 23 in the space of the ipi. In principle, the variational wave function 
needs not to be normalizable at threshold, but this becomes immaterial as 
soon as very long range components are introduced in the expansion flSSp. It 
was checked that the extrapolation method and the direct computation give 
the same result for the frontier. 
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